Construction of spectral functions for medium-mass nuclei 
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This article is aimed at improving the description of lepton-nucleus interactions in the sub-GeV 
energy range. Approximate spectral functions for oxygen, calcium, and argon are constructed and 
used to obtain the electron cross sections in a given scattering angle. Comparison with a sample 
of available experimental data shows satisfactory agreement. Discrepancy between the presented 
model and the systematic computations available for oxygen [O. Benhar et al, Phys. Rev. D 72, 
053005 (2005)] is also found to be very small. Analysis of appropriate kinematical regions leads to 
the conclusion that the obtained argon spectral function should describe well neutrino scattering in 
the 800-MeV energy region. Several approximations used in the model are critically reviewed. All 
the details needed to implement the presented approach in Monte Carlo simulations are given. 
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I. INTRODUCTION 

In the impulse approximation (lA), lepton-nucleus in- 
teraction is described as a two step process: in the first 
step, the lepton interacts with a single bound nucleon, 
and in the second one, the resulting particles propagate 
inside the nucleus. The lA formalism is the basic frame- 
work in which ~l-GeV leptons scattering is described: 
on one hand, it is the approach applied to understand 
electron scattering, and, on the other, it is used to model 
neutrino interactions The status of these two cases 
is quite different: for electrons a lot of experimental data 
exist, whereas for neutrinos precise measurements are 
still missing. Due to this lack of knowledge, reliable theo- 
retical models are needed in the next generation of precise 
neutrino oscillation experiments possibly with liquid 
argon target. 

To construct a successful model of neutrino-nucleus 
scattering, the following procedure seems to be well jus- 
tified: 

(i) relevant kinematical region in energy and momen- 
tum transfer has to be identified, 

(ii) description of a nucleus for electron scattering 
should be formulated in this kinematical region, 

(iii) performance of the electron scattering model must 
be confronted with the existing data, and, if the 
agreement is satisfactory, 

(iv) the same treatment of nuclear effects should be ap- 
plied to neutrino interactions. 

This is the basic logic of this paper, in which we pro- 
pose a model to describe the ~l-GeV neutrino scatter- 
ing off medium-mass nuclei, such as calcium and argon. 
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The article reports continuation of the research started 
in Ref . Q , where a less sophisticated description was ap- 
plied to argon, and Ref. [j], where the model was intro- 
duced. 

In the lA regime, a nucleus is described by means of 
the spectral function (SF). The SF contains information 
about the momentum distribution in conjunction with 
the distribution of binding energy of nucleons inside the 
the nucleus. Evaluation of the SF for medium nuclei 
requires several approximations. In our presentation, we 
try to identify and justify all the theoretical assumptions, 
but the most important argument for the correctness of 
our model is the agreement of its predictions with the 
data for electron scattering. Detailed verification of the 
description is performed using two targets, namely oxy- 
gen and calcium. Oxygen was selected because of addi- 
tional opportunity to compare results with a more sys- 
tematic theoretical approach to modeling of the SF [^], 
whereas the calcium nucleus is most similar to argon's, 
for which precise measurements have been performed. 
Our description of 20 Ca is confronted also with the the- 
oretical results of Butkevich and Mikheyev %]. Finally, 
a comparison with the few known data for electron scat- 
tering off argon Q is done as well. 

Basic computations of quasielastic inclusive scattering 
are standard and can be found elsewhere P. The out- 
come of numerical calculations depends on several as- 
sumptions which specify the implementation of the SF 
model. We apply the recent BBBA05 parameterization 
of the proton and neutron form factors The off-shell 
hadronic current matrix elements are evaluated with the 
use of the standard de Forest prescription [l^. Further- 
more, in the electromagnetic case, we adopt a procedure 
to impose electromagnetic current conservation. Such 
prescription is not unique, and for this reason, in the 
case of weak interactions, we avoid analogous manipula- 
tions with the vector part of the current and rely on the 
de Forest approach only. 

An important ingredient in calculations is the treat- 
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ment of final-state interactions (FSI). There are several 
approaches to deal with them, e.g., the Wentzel-Kramers- 
Brillouin method, relativistic mean-field approach , or 
correlated Glauber approximation ^1^. In the previous 
article 0], we adopted the plane wave impulse approx- 
imation (PWIA) and disregarded FSI beyond the Pauli 
blocking, arguing that this approach is sufficient to de- 
scribe neutrino scattering. In this paper, the presented 
model is validated by confronting it with a large sample of 
electron scattering data and this comparison requires in- 
clusion of FSI. We consider two FSI effects: Pauli block- 
ing and reinteractions of the struck nucleon with the spec- 
tator system described by means of the time-independent 
optical potential [12, [13, Cl- 
in the lepton energy range of '^1 GeV, two dynamical 
mechanisms are most important: quasielastic (QE) scat- 
tering (throughout this article we use the terminology 
of the neutrino community) and single pion production 
through the A excitation. They clearly manifest them- 
selves as two peaks in the electron differential cross sec- 
tion in energy transfer for fixed scattering angle. Our 
numerical computations include the QE dynamics only. 
The reason is that our concern is to define a system- 
atic procedure to construct SF, and it is sufficient to 
test it in the case of the QE process. Thank to this, we 
avoid dynamical issues in which the theoretical situation 
is not completely clear, namely the nonresonant back- 
ground and two-particles-two- holes excitations [13, [13] ■ 
However, we have to pay the price for the constraint on 
the dynamics we adopted; in comparisons of our predic- 
tions with the experimental data, for higher values of the 
energy transfer some strength is systematically missing. 
Strictly speaking, we verify our model only in the kine- 
matical region of energy transfers below the QE peak. 

This paper is organized as follows. In Sec. HH the con- 
struction of our model is described: In Sec. Ill A| we 
present basic formulas for the lepton-nucleus cross sec- 
tion and introduce notation used throughout the article. 
Section III Bl discusses a relation between kinematical re- 
gions in electron and neutrino scattering. In Sec. Ill CI 
a method to approximate SFs is given. The treatment 
of FSI effects in our model is covered in Sec. IIIDI Sec- 
tion III El provides the parameterization of the SFs for 
oxygen, calcium, and argon. The information is detailed 
enough for everybody to be able to implement our results 
in their own numerical codes. 

In Sec, mil our results are compared to large data sam- 
ples for electron scattering off oxygen, calcium, and ar- 
gon, selected according to the conclusions of Sec. IIIBI 
Our predictions are also confronted with other theoreti- 
cal approaches. We observe that the performance of the 
presented description of nuclei is satisfactory and arrive 
at the conclusion that when applied to neutrino scatter- 
ing, the model should produce reliable results. 

Sec.|lV]is devoted to a discussion of the approximations 
used in this paper. We consider plausible modifications 
of the adopted parameters and try to understand how 
uncertain our results for the cross section are. 



Finally, in Sec. |Vl we summarize the conclusions of 
this article. Since the most important features of the 
predictions seem to follow from the very basic assumption 
of the lA, the failure of the model in some kinematical 
situations may be interpreted as a failure of the lA itself. 
The presented results suggest that the lA starts to be 
unreliable when the typical value of momentum transfer 
is smaller than ~350-400 MeV. 



II. DESCRIPTION OF THE MODEL 
A. General information 

We consider quasielastic (QE) electron scattering off 
a (Z, TV) nucleus of mass Ma, which changes its four- 
momentum from k = (i?k,k) to k' = (i5k',k'). Associ- 
ated with this interaction energy and momentum trans- 
fers are uj = E]f^ — E\^i and q = k — k', respectively. When 
the impulse approximation holds, i.e., when only one nu- 
cleon is involved in the primary vertex, nuclear effects 
can be described by means of the spectral function. 

The proton spectral function (SF) of a given nucleus 
P(j3)(p, i?) is the probability distribution of removing 
from this nucleus a proton with momentum p and leaving 
the residual nucleus with energy 

Er = Ma- M + E + TA-1 , 

which includes recoil energy of the residual nucleus 
Ta^i = pV(2Ma-i), compare Refs. [1, [13- The neu- 
tron SF is defined in an analogous way. 

Energy balance of QE production of a free nucleon car- 
rying four- momentum p' = (£'p',p'), 

oj + Ma = Er + Ep' , 

may be rewritten in a usefvil form using removal energy E, 
which is an argument of the SF: 

LU + M-E^TA-i+Ep, . 



In Sec. Ill CI we will justify that the recoil energy can be 
neglected, and therefore from now on, the energy balance 



Lu + M - E = E, 



p' 



(1) 



is used. 

According to the lA, the inclusive electron-nucleus 
cross section is the sum of contributions from protons 
and neutrons: 



da 



dcr (p) ^ da („) 



daj(i|q| dtJiilql d(jj(i|q| 
Each term is expressed by the standard formula 



dLud\q 



El 



EpEpi 



(2) 



:S{Lo + M-E-Ep,)L';:::HZ,t, 
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where the index t denotes the nucleon isospin. The lep- 
tonic tensor is given by 

L^^, — 2{k^k^ + k^ki, — k ■ k g^v) i 

due to neghgible mass of electron, and the hadronic ten- 
sor is 



with the scalar coefficients Hi t and H2.t depending on 



q 



q| and r = —q /{AM ) in the following way: 
Hi,t - r(Fi,t+F2,t)' , 



H2,t — 



The form factors Fi^t = Fi^t{q^) are in turn expressed 
by the appropriate electric Ge, t and magnetic Gm, t form 
factors 

To handle the problem with the off-shell kinematics, 
we use the de Forest prescription treat interacting 
nucleon as free and use free form factors but modify the 
energy conservation to take into account that a part of 
energy transferred by the probe is absorbed by the spec- 
tator system. Comparing Eq. ([T]) to the energy balance 

where the part of energy transfer which goes to the on- 
shell interacting nucleon with p = {Ep,p) is denoted 
by w, one can find momentum-dependent binding energy: 

eB = Ep-M + E . 

Replacing q = (w,q) by g = (w,q) = — eB,q) in the 
hadronic tensor. 



-"cm. t 



-'-'cm, t ; 



(3) 



we obtain the standard description of the off-shell kine- 
matics. 

However, this procedure violates the conservation of 
the electromagnetic current, because q^H^t^-^ t ^ 0. To 
restore it, we have to add a correction to the contraction 
of the tensors: 



with a shorthand notation introduced for 
2k q 



Q = 



V 



|q| , 

141 



When we consider QE muon neutrino scattering, 
its four- momentum is denoted as fc = (_E^,k), four- 
momentum of the produced muon as fc' = (ii'^,k), and 
q= (w, q) = k — k'. The cross section 

da-^ G|cos^.,|q| f,^,s^M^ 



dajd|q| 



E?, 



E„E, 



contains contraction of the leptonic and hadronic tensors 
Lf^i, — '^{kfiki, + k^ki, fc • fc (jf^jy — i^fii/pak^k ), 



where 



Hi 

H2 

H3 



W^4 + 2W + 9>"^)i^5, 



Fj(l+T)+r(Fi+F2)' 



2Fa{Fi+F2), 



1 



= Ho . 
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FaFp - tF'p 



The tensors differ from the ones for electromagnetic inter- 
action due to the axial contribution (in our calculations 
axial mass Ma = 1.03 GeV) and to the fact that, thanks 
to the conserved- vector-current hypothesis, Fi and F2 are 
expressed by differences of the proton and neutron form 
factors; see Ref. j;9|. Considering neutrino interactions, 
we apply the de Forest prescription ([3]) but do not restore 
conservation of the vector current. All the other quan- 
tities are defined and denoted as in the case of electron 
interaction. 



T em tjA"' 
fiiy cm, t 

which is equal to 



TemrrfJ^i^ , remrjM^ 
/XI/ em, t ' ^fiu-^cor, t ' 



,r2 



T cm rrA"' 

^fii/-"- cor, t 



M 



■ClHi 



C2H2. 



The coefficients ci and C2 can be expressed as 
-4(|k||q|Q-L^k-q) 



C2 = ciV"^ + A{uj -Zj)rQ[\<i--p 



p • q 



kq 



(4) 



(5) 



B. Selection of the electron data 

According to the plan outlined in Sec. HI first we iden- 
tify the region in the (w, |q|) plane which is most impor- 
tant for QE neutrino scattering. The energy and momen- 
tum transfers are related to the muon production angle 9 
by the expression 



E^ -uj 



cos W = 



2£;Jk' 



(7) 



where |k'| — ^ {E^, — ujY ~ rnj^. Therefore fixing 9 is 
equivalent to restricting a region in the (oj, |q|) plane. 
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FIG. 1: Analysis of the dependence of the ^^0{u^, fi~) cross 
section on energy- and momentum-transfer. Lines show what 
electron scattering angles 9^ in the process ^®0(e, e') corre- 
spond to the same kinematical region. The standard differ- 
ential cross section in muon scattering angle for ^"0(1^^, /i^) 
is represented by points. 



Points in Fig. [T] show the neutrino differential cross sec- 
tion dcr™°^'^/(i6' for neutrino energy = 0.8 GeV. The 
peak at ^33° is rather broad and ^50% of the cross sec- 
tion comes from 9 S [20°; 56°]. For E^ = 1.2 GeV, the 
maximum moves to ^^22° and the peak becomes narrower 
(not shown in the figure). 

We want to map the allowed kinematical region for 
neutrino scattering, weighted by the cross section, to the 
corresponding region for electron scattering. For electron 
of energy Ee, the relation analogous to Eq. ([7]) reads 



cos^e = 1 



2Ee {Ee 



(8) 



Hence, for a given value of E^ and selected Ee, we can 
map the muon production angle to the electron scattering 
angle: 



To weight the electron scattering angles by the neutrino 
cross section, we calculated the quantity 



d9 da 



weak 



d9e d9 



(9) 



using the oxygen target, described by the Benhar SF. 
(We checked that the Fermi gas model and the effective 
description lead to the same conclusions.) 

From now on, we concentrate on E^, = 0.8 MeV. In 
Fig. [H we show the quantity ([9]) for three selected values 
of electron energy: Ee = 0.88, 1.08, and 1.2 GeV. The 
conclusion is that to describe well the 0.8-GeV neutrino 
scattering, our model should be verified with 1.2-GeV 
electron data at 9e ^ 23°, 1.08-GeV data at Oe ^ 25°, or 
0.88-GeV data at 9e ~ 30°. 

Let us go into more detail. We deduce that for Ee = 
1.2 GeV, the range of muon scattering angle [20°; 56°] 



corresponds to 9e € [15°; 36°] with a maximum at 23°; 
for Ee = 1.08 GeV, to 9e e [17°; 39°] with a maximum at 
25°; whereas for Ee = 0.88 GeV, to 9e e [19°; 50°] with 
a maximum at 30°. The general rule is that for higher 
electron beam energies, the smaller scattering angles be- 
come significant. 

Equation ([5]) is well defined when Ee > E^. For lower 
Ee, this equation may be applied only for the prize of 
a loss of normalization — the form of the denominator ex- 
cludes some of the points in the {lo, ]q]) plane. For ex- 
ample, when Ee — 0.73 GeV is used, 5% of the strength 
is lost and 9 e [20°; 56°] corresponds to 9e £ [22°; 61°] 
with a maximum at 35°. 

In the case of electron scattering off oxygen, the mea- 
surements were performed for scattering angle 32° using 
beam energies 700, 880, 1080, 1200, and 1500 MeV 0,[ll, 
whereas 537- and 730-MeV beams were used for an- 
gle 37.1° [19]. As follows from our analysis, to obtain 
the model which describes well QE neutrino-nucleus scat- 
tering at energy 800 MeV, the most significant electron 
data are those for 880 and 730 MeV. The relevance of the 
experimental points for 1080 and 700 MeV is smaller. 
The energy and momentum transfers which character- 
ize scattering with electron beams of energies 1200 and 
537 MeV are least similar to what is needed but these 
energies are still in the region of interest. The set of data 
for 1500 MeV was collected at too high scattering angle 
for our applications. 

Among a few papers reporting results of electron scat- 
tering experiments with a calcium target [l^, [U, lU 
the most suitable for testing of our model is Ref. 



23|, 



containing data at the lowest scattering angle, namely 
45.5°, in conjunction with the highest values of beam 
energy — up to 841 MeV. We have checked that all the 
measurements at 45.5° correspond to our region of in- 
terest in the (cj, |q|) plane. Obviously, only the data for 
Ee ~ 841 MeV cover the whole region, and the lower 
electron energy is, the more normalization is lost. For 
example, when one uses Ee = 545 MeV, 9 G [20°; 56°] 
corresponds to 9e E [29°; 75°] with a maximum at 46° 
and 27.4% of the strength is lost. Therefore, we rely 
mainly on comparisons with the experimental data for 
higher electron energies. 

Finally we want to explain why we decided to study 
neutrino energy E^, = 0.8 GeV. The reason is that there is 
a lot of relevant electron scattering data to compare with. 
For higher E^, say 1.2 GeV, the situation would be quite 
different — the electron-scattering data at smaller angles 
would be required, but they are missing for the targets 
we are interested in. 



C. How we model spectral function 

The spectral function describes distribution of nucle- 
ons in the (p, E) plane. By integrating out the depen- 
dence on E, the momentum distribution nt(p) is ob- 
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tained: 

nt(p) = j Pt{p,E)dE . (10) 

Our normalization convention is 

J Ptip,E)d^pdE^Nt , (11) 

where the number of nucleons Nt is Z for protons and N 
for neutrons. 

Approximately 80%-90% of nucleons in a nucleus can 
be described as occupying shell-model states and moving 
freely in the mean-field (MF) potential. The rest of them 
take part in interactions. It is natural to decompose the 
SF into the sum of the MF and correlated parts [l^, HB, 

Pt{p,E) = Nt[PriP,E)+Pr"{p,E)] . (12) 

By analogy to Eq. (jlOp . the MF and correlated momen- 
tum distributions are introduced: 

nr(p) = / PriP,E)dE, (13) 

nr' (p) = J Pr"{p,E)dE, (14) 

so the momentum distribution can be described as com- 
posed of two subdistributions: 

n,(p)=iV,[nr(p) +<"■(?)] • (15) 



1. Treatment of the MF part 



widths. Than the MF momentum distribution (fT3l) can 
be expressed as 

-r(P)=E^i'^"(p)i'- (17) 

Let us make the crucial assumption: each level con- 
tributes equally to the MF momentum distribution. It 
means that in Eq. (jl6p for each a we can make the sub- 
stitution 

Co|0o(p)P^nr(p). (18) 
The final form of the MF part of the SF, 

pr(P,i?)=nr(P)^E^"(^"-^) ' (19) 

have to be further specified by the form of the function 
which describes level width. For a given half-width, the 
Breit-Wigner distribution has longer tails then the Gaus- 
sian one, so we found the latter more suitable: 



F^{x) = ^-^e^p{~8xyDl) . (20) 

Therefore we refer to the proposed model as the Gaus- 
sian spectral function (GSF). The factor 8 in the argu- 
ment of exponential function is introduced for further 
convenience. 

Note that the sum in Eq. (fT9|) extends to all occupied 
states. This approach differs from the one presented in 
Ref. and allows to avoid singularities in the argon SF. 

To describe a specific nucleus by its Gaussian SF, one 
needs to know the appropriate MF momentum distribu- 
tion, the values of energy levels, and their widths D^. 



The basic assumption underlying the presented ap- 
proach is the lA, therefore the MF part of the SF can 
be written in the form (compare 0, 0, [1^): 

PI^^{p,E)^Y.^\^-(p)\'MEc+Ta-i-E) , (16) 

with separated contributions from each shell-model 
state a, a ranging from 1 to Nt- Denoting spectroscopic 
factor by Cq,, wave function by (j)a{p), level energy by 
Ea , and a function describing level width by Fa , we have 
omitted the isospin index t for clarity of the notation. If 
interactions between nucleons disappeared, the MF part 
would describe the whole SF (equivalently, all Cq, 's would 
become equal to 1) and each Fa would be the 5 function. 

In this article we are interested in a description of 
medium-sized nuclei, like calcium and argon. Recoil en- 
ergy of the residual nucleus T^-i may then be neglected 
in the MF part of the SF since it is typically ^0.5 MeV 
(see the average MF momenta in Table H]). 

We assume that J Fa{E) dE = 1, what can be physi- 
cally interpreted as the momentum independence of level 



2. Approach to the correlated part 

Interacting nucleons are described by the correlated 
part of the SF. It is a known fact (see Ref. [l^] and ref- 
erences therein) that the two-nucleon interactions domi- 
nate. These short-range correlations give rise to pairs of 
nucleons with high relative momentum. We follow the 
approach of Kulagin and Petti and do not include 
in the considerations interactions of higher order. Than, 
P/^°"' can be expressed analytically in the form: 

prip.E) = <"-(p)o\/^ 

IpI V TT 

X [exp(-apf„i„) - exp(-ap^^^)] . (21) 

The constant a appearing in the above formula is a short- 
hand notation for 3/(4(p^p)/3) with (3 = {A-2)/{A-l) 
and the mean square of the MF momentum (Pmf) de- 
fined as 

(Pmf) = —c my} \ L ' (^^) 
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whereas 



Ir'max 



{/5|P| 



2M(3[E - £;(2) - Ta-i] y 



2M(3[E - £:(2) 



(23) 



The two-nucleon separation energy E'^^'> is an average ex- 
citation of the {A — 2) nucleon system. Since by definition 
averaging should be carried out only over the low-lying 
states, it can be approximated by the mass difference 
£;(2) = Ma-2 + 2M - Ma- 

Because an overwhelming contribution to the corre- 
lated part comes from the peak at 



EKi E 



(2) 



P 

2M 



and the recoil energy Ta~i is less than p2/(2M) by the 
factor {A — 1), therefore Eq. ((23| may be simplified to 



2 

Pmin 



P' 

Jr max 



{/3|P| 
{/3|P| 



2M/3[E - S(2 



2M/3[E - ^(2 



(24) 



For the lightest considered here nucleus, i.e., oxygen this 
simplification yields a < 0.2% change of the cross section. 



D. How we apply FSI 

The struck nucleon moves in nuclear matter and may 
interact with surrounding spectators. Such interactions 
make the nucleon an open system in the sense that mea- 
sured Epi is not equal to its energy in the interaction 
vertex. One can describe this situation in terms of a com- 
plex optical potential, U — V — iW, as proposed origi- 
nally in Ref. jl^. We assume that the potential is time- 
independent. Then the result is equivalent to making in 
Eq. ([2]) the substitution 



S{...) 



W/tt 



W^ + [ V\ 



(25) 



see Refs. [12|, [l3|. The imaginary part of the optical 
potential may be approximated by 



I tj-pl 



(26) 



This article's main interest is a description of medium 
nuclei, such as calcium and argon, therefore the nuclear 
matter density pnuci is assumed to be constant and equal 
to the saturation density psat — 0.16 fm^'^. In the kinc- 
matical region of our interest, the typical proton kinetic 
energy is 100-300 MeV and the nucleon- nucleon cross 
section ajsiN = hiTpp + cTpn) at psat varies between 16.2 
and 19.1 mb [Ig. We set it to the value for 200-MeV 
protons, i.e., to 17.4 mb. 




FIG. 2: Optical potential used in this paper. Dashed line 
represents its imaginary part obtained from Eq. (|26|) and solid 
line the real one from Eq. (|27|l . For comparison, the imaginary 
part calculated also from Eq. ([27} is shown by dotted line. 



The real part of the potential we use is calculated in 
the following way: Reference [s^l gives a Dirac optical 
potential of fpCa fitted to proton-scattering data in the 
energy range 161-1040 MeV as a function of kinetic en- 
ergy of proton and position in the nucleus. Since what 
we need is the potential depending on energy only, aver- 
aging over spatial coordinate should be performed. We 
do it by evaluating the potential at the root mean square 
(rms) radii from Ref. [31[. As a result, we obtain a po- 
tential ?7(p') related to the scalar and vector part of the 
potential in Ref. [s^l by 



i?P' +C/(p') = [M + S{Tp, , rsW + P'' + V{Tp,,fv) • 

(27) 

In the above equation, denotes two parameters, be- 
cause the real and imaginary part of S have different 
values of the rms radius. The same holds true for fy and 
V. For Ip'I > 3.1 fm^^, the real part of U{p') is pos- 
itive, what is inconsistent with the correlated Glauber 
theory [2^. Therefore when |p'| > 3.1 fm~^, we set its 
value to zero, as shown in Fig. [51 

From the few parameterizations of the potential in 
Ref. [lOl we decided to use the one called case 2. We 
checked that the imaginary part of U{p') is then very 
close to W obtained from Eq. ([55)1 (compare the dot- 
ted and dashed lines in Fig. so our approach is self- 
consistent. 

The assumption that the optical potential is time- 
independent leads to folding of the cross section with 
the Lorentzian function [Eq. ([^5)) ] . To cure the resulting 
problem with nonzero cross section for uj < (compare 
Fig. 4 in Ref. [l^), we impose an additional constraint 
on the upper limit of the integration over E, 

E < UJ . 

When FSI effects are not included, this restriction comes 
automatically from the energy-conserving S function. 

As can be seen in Fig. [H the essential effect of the 
imaginary part of the potential is to broaden the QE 
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FIG. 3: Influence of FSI on electron-nucleus cross section. 
Dotted line shows the cross section without FSI, dashed line 
with only imaginary part applied, and solid line with full FSI. 

peak, whereas the real part mainly moves the strength to 
lower oj's. Thanks to these two effects, the agreement of 
the calculated cross sections with the experimental data 
is significantly better. However, the time-independent 
imaginary part of the optical potential overestimates the 
FSf, namely too much strength is redistributed from the 
QE peak to its tails. Wc postpone the discussion of this 
point to Sec. IIVI 

E. Details of implementation 

In this subsection, we want to cover all the details 
of description of three nuclei — oxygen, calcium, and 



argon — by their Gaussian SFs. A procedure to divide 
momentum distributions given in Ref. [32] into the MF 
and correlated parts is presented and justified. We con- 
centrate on n^°"{p) and obtain the MF part from (|15p . 
Then, we show the parameterization of the energy levels 
and comment the way, their widths are obtained. 



1. Momentum distributions 

In Ref. [12, , the total momentum distributions for 
many nuclei are calculated. However, the model de- 
scribed in this paper requires a separation of the MF and 
correlated contributions [see Eqs. (fT9| and ((2T|) ]. Ref- 
erences [13, mi contain plots with momentum distribu- 
tions divided in the way we need. The conclusion from 
these articles is that above |p| = 2 fm^^, the correlated 
part dominates overwhelmingly. We assume that above 
2 fm^^ this contribution is equal to the total distribution, 
and that rif"" (p) may be expressed as the correlated dis- 
tributions given there, i.e., as a sum of two exponential 
functions. Moreover, smooth transition at 2 fm^^ is im- 
posed. 

Distributions from Refs. [13, [13 1 denoted here by 
n{t,p), are calculated up to |p| = 3.585 fm^^. We ex- 
trapolated them smoothly to 5 fm"""^, but it turned out 
to have very little influence on the cross sections. 

The correlated part of the momentum distribution is 
assumed to be of the following functional form: 



(2^r 


Nt 




A 


(27r)i' 


Nt 




A 



for 2.025 fm-i < |p| < 3.585 fm"!, (28) 
for 3.585 fm^i < |p| < 5.0 fm^^. 



In the above equation, A stands for the number of nu- 
cleons. We normalize the momentum distributions intro- 
ducing the factor J-': 



5 fin~ 



47rp2n(i,p) d|p| 



T A 
N't Jo 



To find the values of the parameters in ([281) . we assume 
that ei ^ 62, so that only the e2-containing term is re- 
sponsible for the behavior of nt°"{p) at large momenta. 
By demanding the continuity and smoothness of nl"", 
one can determine C2 and 62 at |p| = 2.025 fm"-'^, while 
C3 and 63 at |p| = 3.585 fm^^. The values of ei are 
taken from [IB]; in Sec. IIVI we will show that ei's do 
not affect the cross sections. The values of Ci are fixed 



by the overall normalization of n™"^, which follows from 
Ref. [12|: the data contained there in Tables II and III 
allow to calculate what fraction of nucleons cannot be 
assigned to any shell-model state and, as a consequence, 
must be described by the correlated part of SF. The nor- 
malization of nf" with respect to nt is the same as the 
normalization of the the correlated SF with respect to 
the total SF. 

Sample outcome of the described procedure is pre- 
sented in Fig. m One can see that at the sewing points, 
IpI = 2.025 fm~^ and 3.585 fm~^, the correlated con- 
tribution is smooth. The total momentum distributions 
and the normalization of the correlated parts alike are 
taken from Ref. ^3^]. Therefore the set of parameters 
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IpI (fm"^) 

FIG. 4: Proton momentum distribution in tElCa from Ref. ^ 
(dots) divided into the MF (dashed line) and correlated part 
(dotted and solid line). Solid line shows extrapolation accord- 
ing to Eq. (|28p with parameters given in Table HI 

TABLE I: Parameters of the correlated part [see Eq. fIE\i ] of 
the momentum distributions from Ref. [32] and their normal- 
ization with respect to the total momentum distribution for 
various nuclei. The last row contains resulting values of the 
average mean field momentum defined in ((22}. 





1o 


40 
20 

proton 


Ca 

neutron 


48 
20 

proton 


Ca 

neutron 




1.0200 


1.0370 


1.0370 


1.0440 


1.0200 


Ci 


2.1280 


4.2150 


4.2700 


4.0040 


4.6700 


ei 


1.4000 


1.7700 


1.7700 


1.7700 


1.7700 


C2 


0.1427 


0.1940 


0.1855 


0.1536 


0.1656 


62 


0.2260 


0.2260 


0.2142 


0.2018 


0.2065 


Cs 


0.1678 


0.2282 


0.2451 


0.2500 


0.2960 


63 


0.2410 


0.2580 


0.2648 


0.2972 


0.2940 


Normal. 


12.00% 


16.20% 


16.20% 


17.10% 


13.64% 


V(p^f) (MeV) 


174.4 


189.1 


187.1 


180.8 


196.4 



collected in Table U can be considered as self-consistent. 

To handle the lack of knowledge of the momentum dis- 
tributions for protons and neutrons in the argon nucleus, 
we apply in the SFs the appropriate distributions calcu- 
lated for onCa. 



2. Description of the energy levels 

In our approach, each shell-model state a is fully 
characterized by two parameters: energy level Ea and 
width Da , defined by means of Eq. ([^0]) . 

The energy levels of calcium shown in Table |ll] result 
from theoretical calculations in Refs. [s^ (for neutrons) 
and [Saffor protons). A few available neutron levels of 
argon [3S] form a pattern very similar to the one of the 
neutron levels of calcium: the distance between lc?3/2 and 
2si/2 is 1.7 MeV for Ar and 1.74 MeV for Ca, whereas 
the distances between 1^3/2 and the Fermi level ap are 
3.5 MeV for Ar and 3.8 MeV for Ca. To reconstruct 
the missing data, we assume that all the neutron levels 



TABLE II: Energy levels 11,133 and widths D„ for toCa. 





protons 

Ea 




neutrons 


Da 


tSl/2 


( .00 


zo 


DD.lZ 


ZO 


lP3/2 


36.52 


15" 


43.80 


15" 


lPl/2 


31.62 


15" 


39.12 


15" 


1^5/2 


14.95 


4" 


22.48 


6" 


2Sl/2 


10.67 


2" 


17.53 


4' 


1^3/2 


8.88 


2'- 


15.79 




ap 


4.71 




12.0 





"Fit to the plots from Refs. [3^.[37t. 

''Our estimate, details in text. 



TABLE III: Same as Table HH but for f^Ar. Details in text. 





protons 

Ea 


Da 


neutrons 

Ea 


Da 


lSl/2 


52" 


25 


62 


25 


lP3/2 


32*" 


15 


40 


15 


lPl/2 


28' 


15 


35 


15 


14/2 


11' 


4 


18 


5 


2Sl/2 


8' 


2 


13.15" 


4 


1^3/2 


6' 


2 


11.45" 


3 


I./7/2 






5.56" 


3 








8.0" 





"Theoretical calculations in 38]. 

''Modified theoretical values for lijCa from [37l| . 



follow the same pattern, see Table UlIl Due to the lack of 
knowledge about the proton levels of argon, we use the 
modified values from calcium. The data for oxygen (39l . 
SO] are collected in Table [IVl 

The widths for most of the calcium levels can be de- 
termined by fitting to the plots of energy distribution in 
papers [H, [13]. We estimate the remaining ones using 
the fact that Da should be, approximately, a function of 
a distance from the Fermi level ]36, 3Zj JiSj 41,] : 



" {Ea-Epr+a^' 

To get ZJq's for argon, we assume that their values lie on 
roughly the same curve as for Ca. 

We did not find the energy distribution for oxygen, cal- 
culated in the same way as for calcium in Refs. [36, 37|. 
Since oxygen nucleus plays only the role of a testing 



TABLE IV: Same as Table HT] but for ^|0. The values of D, 
are obtained differently; see details in text. 





protons 

Ea 


Da 


neutrons 

Ea 


Da 


lSl/2 


45.00 


70 


47.00 


70 


lP3/2 


18.44 


4 


21.80 


4 


lPl/2 


12.11 


4 


15.65 


4 
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FIG. 5: Cross sections of the process ^^0(e,e') at miscellaneous values of beam energy for scattering angles 32° 0,13 ''■nd 
37.1° [l^. Results for the GSF (solid line) are compared to the Benhar SF 5] with the same FSI (dashed line) and the Fermi 
gas model without FSI (dotted line). The values of |q| at the peaks are 637 MeV (for beam energy 1200 MeV), 573 MeV (for 
1080 MeV), 466 MeV (for 880 MeV), 371 MeV (for 700 MeV), 441 MeV (for 730 MeV), and 325 MeV (for 537 MeV). 



ground for our model, we decided to obtain the proton 
DaS directly from the energy distribution in the Ben- 
har SF, and use the same values for neutrons. Thus we 
avoided additional discrepancies between the two descrip- 
tions. 



formula 



|q| = ^ 2E]^{E]^ - uj){\ - cost 



This value depends rather weakly on uj and therefore it 
provides quite good characteristics of the whole peak. 



III. RESULTS 

A. Electron scattering 

The goal of this subsection is to confront the model 
presented in Sec. |TT] with the existing electron scatter- 
ing data. Since description of the dip region and the 
A excitation is ambiguous [1, H^, , our considerations 
include QE interactions only and we test predictions of 
the obtained SFs in energy transfers below the QE peak. 
Figures [5H7| show comparison with a wide spectrum of 
experimental points. The missing cross section for en- 
ergy transfer above the QE peak may be attributed to 
the two-nucleon interactions, A production, and nonres- 
onant background. In the captions, we give the momen- 
tum transfer at the QE peak calculated according to the 



1 . Oxygen 

We start with the oxygen target. Figure [5] presents 
the predictions of three models and the data from 
Refs. 0, E m. The dotted line corresponds to the 
Fermi gas (EG) model (Fermi momentum = 225 MeV, 
binding energy es = 25 MeV, no FSI), the solid line 
shows the cross sections of the oxygen GSF with FSI 
as in Sec. IIIDl and the dashed line depicts results for 
the Benhar SF with the same FSI. Differences between 
our model and the more systematic SF are of the size 
of the error bars. The main source of these differences 
is another momentum distribution, see Fig. [TT] and the 
discussion in Sec. IIVI Both SFs reproduce the shape and 
height of the QE peak quite well, but underestimate the 
cross section at low w's. This discrepancy may be at- 
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FIG. 6: Cross sections of ^''Ca(e, e') scattering at angle 45.5° and miscellaneous values of electron beam energy [2^ . Calculations 
for the GSF (solid line) are compared to the results of Butkevich and Mikheyev (dashed line), and the Fermi gas model 
(dotted line). The corresponding values of |q| at the peaks are 602 MeV (for beam energy 841 MeV), 561 MeV (for 782 MeV), 
531 MeV (for 739 MeV), 490 MeV (for 681 MeV), 453 MeV (for 628 MeV), 395 MeV (for 545 MeV), 342 MeV (for 471 MeV), 
297 MeV (for 408 MeV), and 254 MeV (for 350 MeV). 



tributed to the unsatisfactory treatment of FSI effects, 
because they tend to increase the cross section in this 
region. The best agreement with the data is obtained for 
the 880-MeV electron beam, whereas the worst one cor- 
responds to i?k — 537 MeV. Fortunately, the latter set of 
data is least relevant in our analysis (see Sec. Ill Bj ). 



2. Calcium 

For calcium, we compare in Fig. [5] the cross sec- 
tions obtained using the FG model {pp — 249 MeV, 
es = 33 MeV [13]; represented by the dotted hue), the 
GSF (solid line), and the calculations of Butkevich and 
Mikheyev §\ (dashed line) to the sample of electron scat- 
tering data collected at scattering angle 45.5° and various 
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FIG. 7: Left panel: Comparison of the cross section of GSF (solid line] and the FG model (dotted line) with experimental 
points for Ar(e,e') at beam energy 700 MeV and scattering angle 32° 7]. Right panel: Same, but for oxygen. Note that in 
both cases the similar accuracy is obtained. The value of momentum transfer at the peaks is 371 MeV. 



beam energies (23[. Only our model includes FSI effects. 

The FG model describes very well the position and size 
of the QE peak for the highest values of beam energy. 
When the energy is lower than 700 MeV, it obviously 
fails. 

Despite the fact that the approach of Butkevich and 
Mikheyev @] is based on the SF, it yields results very 
similar to these for the FG model when oj is near the value 
corresponding to the QE peak or higher. The reason of 
this behavior lies in too simple treatment of the MF part 
in their SF: the energy distribution was limited to a single 
5 function. 

For energies 628-841 MeV, the accuracy of the GSF is 
very good. The occurring discrepancies can be explained 
as a contribution from the A production. At small values 
of energy transfer, the cross section is slightly overesti- 
mated. It means also that the QE peak is slightly under- 
estimated, because FSIs based on a folding function do 
not change the total cross section. Note that the agree- 
ment with the data in our region of interest (see Sec. Ill B| ) 
is better than in the case of oxygen. It may be attributed 
to the way FSI is introduced: Density of nucleus is as- 
sumed to be constant and equal to the saturation density 
of nuclear matter; this approximation should work better 
for heavier nuclei. The real part of the optical potential 
used in our computations was obtained for calcium and 
should work better for this target than for oxygen. 

For i?k < 545 MeV, our model fails to describe the 
position and shape of the QE peak. However, the inac- 
curacy of the FG and the approach of Ref. is visibly 
more severe. Similar problems for oxygen occur when 
< 700 MeV. At the first glance, there is no connec- 
tion between these two cases. But when we have a closer 
look at the values of the momentum transfer at the QE 
peak, we will discover that our model starts to lose ac- 
curacy when momentum transfer is lower than ^^350- 
400 MeV. It could be related to simplifying assumptions 
of our approach, treatment of FSI effects or the very basic 
assumption — the lA. The models [5, i23] based on differ- 
ent from our approximations apart from the lA and with 



more systematical treatment of FSI suffer a similar draw- 
back. It suggests that it is the loss of reliability of the 
lA, what is responsible for the problem. 



3. Argon 

Before applying our model to neutrino interactions, we 
perform the final test by confronting it with the data for 
electron scattering off argon. We have found only one 
such experiment [7[ , which measured the cross section of 
700-MeV electrons scattered at 32°. 

In the left panel of Fig.[7l predictions of the argon GSF 
and the FG model [pp = 251 MeV and es = 28 MeV) 
are presented. The accuracy of the GSF is clearly better 
than that of the FG model. The result for our model, 
shown by the solid line, does not describe properly only 
the cross section at very low values of energy transfer. We 
have faced the same problems for oxygen and calcium 
and interpreted it as a breakdown of the lA at |q| < 
350-400 MeV. In the considered case of scattering off 
argon, the momentum transfer at the QE peak is equal to 
371 MeV. When we compare the result for argon with the 
one for oxygen in exactly the same kinematical conditions 
(see right panel of Fig. [7]), we can see that the level of 
accuracy is comparable. The same holds true also for 
comparison with scattering off calcium for electron-beam 
energy 471 or 545 MeV. Therefore, we expect that if the 
typical |q| was higher, the agreement with the data for 
argon would be better. 

We have observed that even for argon, the neutron 
SF may be approximated by the corresponding proton 
SF, as far as electron scattering is concerned. It can be 
explained by the fact that the contribution of neutrons to 
the inclusive cross section is small, what suppresses the 
differences between the SFs. This contribution is equal 
to 13% for 700-MeV electrons scattered at 32° and rises 
to 23% when the beam energy is increased to 1200 MeV. 
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E,, (McV) 

FIG. 8: Quasielastic differential cross section da^'^'^/dE^, of 
fgAr as a function of produced muon energy E^^ for the GSF 
(solid line), approach of Ref. [31 (dashed line), and the FG 
model (dotted line). 



B. Neutrino scattering 

In the case of neutrino scattering, quantities of inter- 
est are the total cross section and the differential cross 
section in = —q^ or in energy transfer (equivalently: 
in energy of produced muon) . 

Figure [8] depicts differences between da^'^'^^ / dE ^ for 
the argon GSF (solid line) , the SF we described in Ref. Q 
(dashed line), and the FG model (dotted line). One can 
see that the SFs introduce significant reduction of the 
cross section, mainly in the region of low energy trans- 
fers. The line representing the predictions of the GSF 
model is slightly wiggly, because when uj increases, lower- 
lying energy levels consecutively start contributing to the 
cross section. There are no singularities in the cross sec- 
tion, and in this sense, the GSF is more realistic then 
the SF from Ref. [3|. Effects of FSI are not taken into 
account except those from Pauli blocking, but their in- 
fluence on the cross section da^'^'^^/dE^ is rather small 
(see Fig. 14 in Ref. 0] showing the impact of introducing 
FSI on da^'^^^/dQ^ ). The purpose of Fig. [5] is to show 
discrepancy of our description of argon nucleus and the 
FG model, commonly used in Monte Carlo simulations. 

The results for neutrinos cannot be directly confronted 
with experimental data. Therefore, we first identified, in 
Sec. IIIBl the region in the (w, |q|) plane which is most 
important for the 800-MeV neutrino scattering. Than 
we substantiated accuracy of our approach: we showed 
in Sec. IIII Al that it describes well kinematical aspects of 
nuclear effects. This whole analysis allows us to expect 
that using the presented approximation of the SF, we 
model neutrino interactions at a similar level of accuracy 
as achieved in the case of electron scattering. 



IV. DISCUSSION OF PRECISION 

Our approach is based on many approximations and in 
this section, we would like to understand how uncertain 




FIG. 9: Dependence of the cross section on the form fac- 
tors. The dipole parametrization (dashed line) produces ~3% 
higher result than the BBBA05 one [^. 



|[jCa(e,e'), 628 MeV, 45.5° 
40 I 




uj (MeV) 

FIG. 10: Influence of the procedure restoring current con- 
servation [Eq. (|4])] on the cross section. The result obtained 
without it is represented by the dashed line. 



our final predictions are. 



1. General remarks 

a. Form factors. Different choices of parameteriza- 
tion of the electromagnetic form factors may change the 
results by a few percentages. As shown in Fig. [9l the 
dipole parameterization yields the cross sections higher 
than the BBBA05 one Q used in this article. The dis- 
crepancy at the QE peak is ~3.2% for beam energy 
350 MeV and ~3.5% for 841 MeV. 

b. Current conservation. Describing both electron 
and neutrino interactions, we applied the de Forest pre- 
scription to describe the off-shell kinematics. However, 
this leads to a loss of conservation of the electromagnetic 
current in electron scattering and of the vector current 
in neutrino interactions. The procedure by which we 
restore it in the electron case, modifies the cross section 
mainly above the QE peak; see Fig.fTUl For beam energy 
628 MeV, the effect is as small as 1.5% at the peak and 
it decreases when the energy becomes larger. 
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FIG. 11: Intrinsic inaccuracy of our model arising from the 
treatment of the MF part of the SF. Calculation for the Ben- 
har's exact SF of oxygen (dashed line) are compared with 
result for the GSF with the same momentum distribution. 



FIG. 12: Uncertainty of the the cross section with respect 
to used momentum distribution. Solid line shows result for 
momentum distribution from Ref. [3^ fused throughout this 
paper) and the dashed one from Ref. [Sa ]. 



c. Simplifications in the mean-field SF. In the 
derivation of Eq. (I19p , we made two simplifying assump- 
tions: level widths do not depend on momentum and 
contribution to the momentum distribution of each level 
is the same [see Eq. p^ ]. Figure [TT] illustrates the loss 
of accuracy due to these simplifications. To depict their 
influence, we use the momentum distribution calculated 
from the Benhar SF instead of the one from Ref. [3^ . 
Since the level widths of oxygen are obtained by fitting 
to the energy distribution of the Benhar SF, a slightly 
different shape of the predicted QE peak is the result 
of the simplifying assumptions only. We checked that 
for other values of beam energy discrepancy does not 
increase. Therefore, we conclude that the GSF can be 
considered as quite good approximation of the more sys- 
tematic approach. 

d. Parameterization of the momentum distribu- 
tions. Application of the momentum distributions from 
Ref. |3^] in our model requires dividing each of them into 
the MF and correlated parts. It involved introduction of 
a few parameters. To find out how much choice of these 
parameters influences the cross sections, we calculated 
first n^:"" for the oxygen normalized as in Table [H but 
with ei = 1.770 (instead of 1.400): 

<"(P) = ^y[2.670exp(-1.770p2) 

+ 0.2128 exp(-0.303 p^)] 

for < IpI < 2.025 fm^^. When it is apphed, the cross 
sections change less than 0.1% in the considered energy 
range. Thus, we do not need to pay much attention to 
parameter ei, as far as the same normalization is kept. 
Second, we found the distribution with ei ~ 1.400, but 
with the normalization 16.2% (instead of 12.0%): 

<"(P) = ^y[3.9228exp(~1.400p2) 

+ 0.0736 exp(-0.091p2)] 




E,, (McV) 

FIG. 13: Same as Fig. [T^] but for scattering off argon. 



at the interval < |p| < 2.025 fm^^. The above distri- 
bution leads to the cross sections changed by up to 2.2%. 
We have analyzed a few such modifications and in each 
case we have found that the influence of the normaliza- 
tion is greater than that of ei. It is because variation of 
the parameter ei only redistributes the strength within 
given part of the momentum distribution [and as a con- 
sequence modifies the parameter a in Eq. (j2ip ]. whereas 
variation of the normalization changes the way some part 
of the strength is treated. 

e. Momentum distributions. Both for oxygen and 
calcium, the momentum distributions are given by ana- 
lytical formulas in Ref. [35^1 . In Fig.[T^ we show that even 
though they predict slightly higher QE peak, the yielded 
cross section is lower. Because the calcium momentum 
distributions are used for argon, its description "inherits" 
the same uncertainties, see Fig. [T31 Throughout this pa- 
per, we rely on the distributions from Ref. [32], because 
they are obtained in more systematic calculations. 
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FIG. 14: Estimation of the uncertainty due to unknown mo- 
mentum distribution of argon. The cross section calculated 
using the SF of igAr with momentum distribution of ItJCa 
(solid line) and loCa (dashed line). 
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FIG. 15: Influence of the level width on the cross section. 
Solid line: calculation with the values from Table Hill Dashed 
line: values multiplied by 3. 



2. Case of argon 

In addition to the already described sources of uncer- 
tainty, the description of argon nucleus suffers from the 
lack of the available momentum distributions and knowl- 
edge of energy levels. We estimate them using the infor- 
mation for loCa. For this reason, a few words of comment 
on the accuracy for this specific nucleus are needed. 

a. Momentum distributions. The surplus neutrons 
modify both the proton and neutron momentum distri- 
butions. A similar situation appears for |§Ca, where the 
distributions are available [33|. We have used the '^G& 
momentum distributions to estimate how these modifi- 
cations can affect the argon cross sections; see Fig. [TH 
The proton cross section was increased by 4% and the 
neutron one was decreased by 3.8%. The overall increase 
is equal to 2.9%. The number of surplus neutrons in fgAr 
is smaller than in IgCa, therefore we expect this effect to 
be smaller too. 

h. Level widths. Due to the lack of any knowledge 
about the level widths of argon, we use the values for 




FIG. 16: Comparison of the cross sections obtained with 
the energy levels from Table HITl (solid line), levels shifted by 
+5 MeV (dotted line), and by -5 MeV (dashed line). 

calcium. Figure [T51 presents that DaS three times larger 
than those given in Table IIIII change the cross sections 
only up to 2% (decrease at the peak). Narrower levels 
gives barely noticeable difference: 0.23% for the widths 
divided by 3, and 0.53% for divided by 100 (increase at 
the peak). For da^'^'^^ / dE^, the more the levels overlap, 
the less wiggly the cross section is. 

c. Energy levels. The argon energy levels may differ 
from the used ones. We may expect that the discrepan- 
cies in Table IIIII are distributed randomly, and so a part 
of their influence on the cross section is diminished. Fig- 
ure[Tn]shows that even if every level is shifted by the same 
value, chosen to be 5 MeV, the cross section does not 
change dramatically — the QE peak only moves a little 
bit. We conclude that the way to increase the accuracy 
of the presented argon SF is to apply the actual values of 
the energy levels; the degree in which they are smeared 
has minor influence on the cross section, especially in the 
case of electrons. 



3. Final state interactions 

a. Real potential. To find out if one can approxi- 
mate the real part of the potential by a constant, we 
have applied the value 10 MeV. In the case of oxygen, it 
slightly improved agreement with the experimental data. 
However, the same value employed to calcium decreased 
the level of accuracy of the model. It might suggest that 
the real potential for oxygen is deeper than the potential 
shown in Fig. [51 

b. Imaginary potential. The use of the imaginary 
part of the potential ?7(p') defined in Eq. (P7)) instead of 
approximation (|26p has minor influence on the obtained 
cross sections. Typical change is a ^1%-increase. We 
conclude that for practical purposes these two approaches 
are equivalent. 

c. Cross section. When evaluating imaginary po- 
tential , we have fixed the nucleon-nucleon cross sec- 
tion to 17.4 mb, which corresponds to nucleon kinetic en- 
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ergy 200 MeV. In principle, one should take into account 
the cross section's dependence on energy. Therefore, to 
check validity of our approximation, we have used the 
exact nucleon-nucleon cross section in the energy 
range 100-300 MeV, most important for the discussed 
kinematical region. For 545-MeV electron scattering off 
calcium, the result decreases by 1.1%. When beam en- 
ergy is higher, the effect is even smaller. 

d. Density of nucleus. We have assumed that the 
density of nucleus is equal to the saturation density, de- 
spite the fact that in reality its average value is smaller. 
However, the quantity of interest is not the density it- 
self but rather PnuciCAw. This product decreases by 7% 
(15.4%) when pnuci changes to 0.14 fm^^ (0.12 fm^^), 
i.e., by 12.5% (25%). Since the corresponding increase of 
the electron cross section is only 1% (2.4%) at the QE 
peak, our approach seems to be well justified. 

e. Folding function. Employing Lorentzian folding 
function, i.e., neglecting correlations between nuclcons in 
nucleus is a crude approximation |l2 ^ 3, . Compari- 
son to the results presented in Ref. 5| suggests that an 
accurate approach could yield the cross sections higher 
at the peak by up to ~15% and with lower tails. Pre- 
cise comparison is difficult because in our computations 
contribution of the A resonance is missing. 



the (w, |q|) plane which is most important for neutrino 
quasielastic interaction. The presented model to describe 
nuclear effects has been then tested using the electron 
scattering data which lie in this region. The obtained 
agreement is good in the case of oxygen and very good 
for calcium. Moreover, our approximation reproduces 
results of the Benhar SF for oxygen with a satisfactory 
degree of accuracy. Detailed discussion of uncertainties 
due to many simplifications of our model have lead us to 
the conclusion that all of them are of the order of a few 
percentages. 

In addition, we have observed that when the typi- 
cal value of the momentum transfer is less than ^--^350- 
400 MeV, systematic discrepancies between the presented 
model and the electron data occur: the shape of the cal- 



culated cross section 



is not suitable to fit the data 



and also increasing amount of strength is missing at low 
energy transfers. A similar problem is present in other 
models 0, [281] what suggests that the source of the prob- 
lem may be the loss of reliability of the impulse approx- 
imation. 

In this paper, we tried to give all the ingredients used 
in our numerical computations to allow implementation 
of our spectral functions in neutrino Monte Carlo gener- 
ators. 



V. SUMMARY 

The main goal of the article is to improve description 
of neutrino scattering off argon in the 1-GeV energy re- 
gion. We have presented the way to calculate approxi- 
mate spectral functions of medium nuclei and applied it 
to electron scattering off oxygen, calcium, and argon tar- 
gets. For neutrino interactions precise experimental data 
are missing. Therefore, we have identified the region of 
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